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Abstract. The pion multiplicity distribution is widely believed to reflect the statistical 
aspects of pp annihilation at rest. We try to reproduce it in a grand canonical picture with 
explicit conservation of electric charge, isospin, total angular momentum, and the parity 
quantum numbers P, C, and G via the projection operator formalism. Bose statistics is 
found to be non-negligible, particularly in fixing the interaction volume. The calculated 
pion multiplicity distribution for (n n ) = 5 turns out to depend strongly on the conserva- 
tion of the angular momentum and connected quantum numbers, as well as on the spin 
state occupation in S-wave annihilation. However, the empirical Gaussian pion multiplic- 
ity distribution cannot be reproduced. This calls in question either the statistical ansatz 
or the rather old data themselves. 
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1 Introduction 

Antinucleon-nucleon annihilation is dominated by nonperturbative QCD effects, and up 
to now there exists no consistent description of its microscopic dynamics within QCD due 
to the complexity of the underlying mechanisms. It turns out, however, that a number of 
global features can be reproduced at least qualitatively by phenomenological approaches 
(for a review, see [[I], Q). In particular, statistical models in several variations have been 
applied to the pp annihilation system (for a review of the early works, see 0), and even 
the fireball picture has been used motivated by the Maxwell-Boltzmann like form of the 
pion-energy spectra 

Especially the pion multiplicity distribution has always been interpreted as a strong 
hint for a statistical behaviour of the pp annihilation system: the data 0] seem to be fitted 
by a Gaussian distribution 



where A = 14 and a = 1.6. It is therefore natural to try to reproduce the empirical pion 
distribution by using a statistical model. 

On the other hand, conservation of the characteristic quantum numbers of the pp 
system, like isospin [0]-|§ or total angular momentum ||10||-||13||, are known to play an im- 
portant role and to severely constrain the statistical phase space. In this paper we address 
the question how these conservation laws influence the pion multiplicity distribution in 
pp annihilation at rest. Comparison of our calculations with the data is hampered by a 
serious experimental problem: the occupation of the various spin-states J PC of the proto- 
nium atom just before its annihilation is apparently not distributed in a statistical way, 
but depends on the density of the hydrogen target and cannot be experimentally con- 
trolled up to now fT3| , |T5p . S-wave annihilation occurs from either a spin singlet (J = 0) 
or a spin triplet (J = 1) state, but with unknown relative occupation. Therefore, we here 
study systematically the gross features of S-wave pp annihilation (e.g. the pion multiplicity 
distribution) as a function of the ratio 




(1) 



with D = 0.9 and (n n ) = 5. Fig. [T| also shows in addition a modified Poissonian distribu- 
tion, which fits the data equally well, 




(2) 
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2 THE MODEL 




Figure 1: Pion distribution of 
pp annihilation at rest. The em- 
pirical data (histogram) [4] are 
compared with a Gaussian distri- 
bution (full line) and a modihed 
Poissonian (dashed line). 



between the spin singlet and triplet occupation probabilities. 

In our approach, we will use a statistical model in connection with the projection 
operator formalism which allows for the explicit conservation of isospin as well as of the 
so-called external quantum numbers, i.e. the angular momentum and P, C, and G parities. 
Instead of the usual Boltzmann approximation, Bose statistics is employed (for reasons 
that will be discussed) and all mesons and resonances up to the ai(1260) are included 
in our calculations. At a fixed average pion multiplicity of (n,) = 5 we discuss the de- 
pendence of the volume parameter on temperature and then present the pion multiplicity 
distributions for several values of X r . 



2 The Model 

It is rather easy to implement baryon number and strangeness conservation into the 
statistical approach: one builds the total partition function from multi-mesonic states 
only, each of which has vanishing strangeness S = 0. Conservation of the electric charge 
Q is included automatically once the third component of the isospin is conserved, i.e. 
I3 = 0, and the constraining effects of P, C, and G parity conservation can easily be taken 



into account |L3| if the angluar momentum is fixed, which together with the total spin 
is subject to the conservation of the total angular momentum of the system. So, it will 
be sufficient to set up a procedure which allows for the projection of the total partition 
function onto the isospin quantum numbers (J, I3) and the angular momentum L. 
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A convenient and, above all, consistent way of treating conservation of both Abelian 
and non- Abelian quantum numbers in the context of a statistical model is provided by the 
projection operator formalism. The simplest procedure would be to calculate the grand 
canonical partition function Z gc of the produced mesons, which are supposed to behave 
like free particles in a spherical "box" while the interaction between these mesons is 
accounted for by resonance production |T6| , 17 1. The unphysical multi-meson states which 



are contained in Z gc due to quantum number fluctuations then have to be projected out by 
suitable projection operators. This is a well established method JT^| - |21] which has been 



mostly applied to internal quantum numbers like £7(l)-baryon number or -strangeness and 
SU (2)-isospin, but can also be extended to treat "external" quantum number conservation 
like that of angular momentum fl3 |. Following this line, the trick is not to calculate Z gc 
itself but a so-called generating function Z. In the case of isospin and angular momentum 
conservation (AMC), the generating function can be written as 

Z(T, V, di, di) = tr [exp(— j3H + idi ■ I + idi ■ L)\ , (3) 

with the symmetry group parameters dj and d^ for the SU (2)-isospin and S'0(3)-rotation 
group, respectively. It can be easily shown, that this generating function can be expressed 
in an alternative way, namely 

z{t, v, d h d L ) = J2 { J^ {2 1\ 1} • ■ x L (d L ) . (4) 



Here, Xi and Xl are the group characters [22], and the summation has to be done over 



all irreducible representations of the symmetry groups under consideration. Zj l is a con- 
strained partition function which contains only those multi-meson states which are eigen- 
states of I 2 and L 2 with eigenvalues /(/ + 1) and L(L + 1) and which is the actual 
quantity in demand. In order to project Z on Zj t L, the orthogonality relation of the group 
characters has to be wrapped up into a corresponding projection operator. 

Let us now write down some explicit expressions for the general case of r different 
particle species being indexed by j = 1, . . . , r. In particular, the particle species are 
distinguished by their isospin t, third component £3 and mass, i.e. the index j is only 
a short notation for (t,t^,m). Furthermore, two particles are considered to be identical 
if they correspond in j, I and their one particle energy e(k). Because the geometry of 
the quantum gas resulting from pp annihilation at rest is supposed to be approximately 
spherically symmetric, the projection onto angular momentum eigenstates should not 
prefer a certain direction. Hence, no projection on the third component of the angular 



momentum is needed, and for this reason we are allowed to replace L by L 3 [22] in 
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2 THE MODEL 



In order to get the generating function of the system, the trace in equation (^) is worked 
out with eigenstates of (i 2 , 1 3 ) and (L 2 , L 3 ), and the occupation number representation is 
employed (for details, see 



Z(T,V,a,u>) = tr [exp(— (3H + idl + iuoL 3 



r oo oo 

nnn 

3=1 1=0 k t =0 



(^i 4 («)) nUi -(^H) n " i - ex pK^-K)) 



(5) 



where (3 = 1/T and qL is the relativistic energy of a particle of type j. The occupation 
numbers n\ fc; denote the number of particles of type j in a state with angular momentum / 
and magnitude of momentum k\. Note, that the momentum k depends on I and the inter- 
action volume V, because the radial component of the angular momentum eigenfunction 
is a spherical bessel function, ji(kr), which must vanish at the edge of the spherical "box" . 

Usually, at this stage the Boltzmann approximation is invoked. The argument is || 
that Bose statistics makes at most only a ten percent correction in the case of the one 
pion partition function, and the correction decreases rapidly with increasing particle mass. 
However, we will show that when going to multi-mesonic states then the error resulting 
from the Boltzmann approximation can be much larger and Bose statistics must be taken 
into account. Thus, after some simple transformations and proceeding with Bose statistics, 
equation (H) can be rewritten in the following form: 



Z(T, V, a, uj) = exp 



oo r i 



la 



n=l j=\ 



1=0 



(6) 



Due to Bose statistics, more than one particle of the same type can occupy a single level 
characterized by its angular momentum eigenvalue I and its energy e{k{). In the above 
expression, the partition function for n particles of type j, each with the same angular 
momentum I, is given by 



h=o 



(7) 



where the sum is over all energy eigenvalues for the given value of I (obtained as zeroes 
of the spherical Bessel function, see above). The classical Boltzmann approximation is 
contained in (|6]) by considering only the term n = 1. 

In general, the form of the projection operators depends on the choice of the symmetry 
group parameters di and a^. For convenience we choose for S"f/(2)-isospin the Euler angles 
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(a, ft, 7) and for SO (3)-angular momentum the unit vector n of the rotation axis together 
with the rotation angle u. Then, we can write for the isospin projector 

A,7 8 =o = Jdajd(3jd 1 sin/3P 7 (cos/3) , (8) 

000 

where the group character Xi was replaced by the matrix element £)g (a,/3, 7) = 
P/(cos(/3), because we want to project onto (I, ^3 = 0) rather than only /. In the case of 
AMC we get 

7T 

or ii« p 

P L = -l^T j dnj ^sin(f)sin((2L + l)f) . (9) 


The only way to apply these projectors directly to Z of equation (Q) is to expand the 
exponential and the resulting expressions until we have a sum over all possible final multi- 
meson states contained in the generating function. Because of the huge amount of such 
states (all mesons up to the ai(1260) are included in the calculations), the evaluation 
is done by the computer. In particular the partition functions z J nl (T,V) are evaluated 
numerically. Furthermore, the evaluation remains treatable only up to a total number of 
N = 7 mesons. This cut-off in the partition function can be well justified by looking at 
the data: the contribution of channels with N > 7 is negligible. 



3 Average Pion Mulitiplicity and Multiplicity Distri- 
bution 

In a first step, we want to adjust the model parameters T and V to the empirical pion 
expectation value (n w ) = 5. In addition to the multiplicity of the directly produced pions, 
one must also calculate the multiplicities of all the resonances in the system which then 



decay according to [23] into "secondary" pions. 

In order to extract multiplicities from the total partition function of the system 
Z(T,V, B, S, I, . . .), we introduce for every particle type j a chemical potential /ij, and 
every partition function z J n t has to be multiplied by the factor 

(X j ) n = expft^-n . (10) 

Then, the particle mulitplicities result from the derivative of the logarithm of the total 
partition function with respect to the fugacities at X J — 1, 



(N J ) = d x An Z(T, V,B, S,I, . . . , X j , . . .) 



(11) 
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3 AVERAGE PION MULITIPLICITY AND MULTIPLICITY DISTRIB UTION 



In Fig. ^| the interaction volume of the pp annihilation system is plotted against the 
temperature of the system for a fixed value of (n n ) = 5. The dashed lines correspond to 
(/, I^) conservation only. Here, the deviation of Boltzmann statistics in curve (a) from 
Bose statistics (b) is not very significant. However, if AMC and the conservation of the 
parity quantum numbers are included in our calculations, then the difference is striking: 
Bose statistics (d) allows for a rather small and therefore much more realistic volume in 
a temperature range of T = [140 . . . 200] MeV whereas the Boltzmann case (c) yields a 
reasonable volume only for very high temperatures of T > 200 MeV. This is a significant 
phenomenological improvement compared to previous studies which were plagued by un- 
realistically large fireball volumes ^ [|. This effect is obviously an intricate combination 




Figure 2: Dependence of the interaction 
volume V on temperature T for (n n ) = 
5 for various scenarios: (a) Boltzmann 
statistics, isospin conservation only; 
(b) Bose statistics, isospin conserva- 
tion only; (c) Boltzmann statistics, ad- 
ditional AMC and conservation of C, 
P , and G parity; (d) Bose statistics, 
additional AMC and conservation of C , 
P , and G parity. 



200 180 160 140 120 100 
T [MeV] 



of both Bose statistics and conservation of external quantum numbers: the Boltzmann 
curves (a) and (c) show opposite behaviour, compared to the Bose curves (b) and (d), 
when going from the dashed to the full lines, i.e when including AMC and conservation 
of the parity quantum numbers. Note, that in (c) and (d) we did not distinguish between 
the protonium spin states, because in our calculations the volume is nearly the same (up 
to 5 percent) for _+ and 1 . This means, that in our model the volume is independent 
of the spin-state ratio X r . 

Once the parameters T and V are known, the pion multiplicity distribution can be 
calculated in a second step. To this end we use the following trick: After expanding the 
projected partition function into a sum over all the allowed multi-meson channels, we 
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£ PM , (12) 



multiply each channel by 



where pin^) is the probability (calculated from the data tables |23]]) that, after all reso- 
nances have decayed, this particular channel yields final pions. Then, in order to be 
able to project out a particular final state with n n pions, we multiply each term in the 
sum on the right hand side of eq. flT2| ) with a Lagrange multiplier f(n n ). Explicitly, we 
make the replacement 



oo 

Z n u h ' ' ' Z n r ,l r 1 > 22 Z n lt h ' ' ' Z n r ,l r ' Pni,..'.% r ( n n) " f ( n n) , (13) 
n^=0 



where jPn{"^n ( n ^) * s ^he probability for having n n final pions produced by a channel 
containing n\ resonances of type ji, resonances of type j'2, etc. If the logarithm of 
the so prepared partition function is derived with respect to f(n n ) (setting f(n n ) = 1 
afterwards), we get the total probability of having pions in the final state as 

PM = d fM In \Z(T, V,B,S,I,..., f(n n ), . . .)} . (14) 

The calculated pion distribution for the simplest case of isospin conservation only in the 
Boltzmann approximation is shown in Fig. for two extreme temperatures, T = 100 MeV 
and T = 200 MeV. The main effect of high temperatures is to give the high mass 
resonances more weight in the partition function, and the resulting pion distribution 
should be broader than at low temperatures. Especially for Boltzmann statistics, however, 
this effect seems to be rather small, see Figure ^|. The Bose case is a little bit more sensitive 
to the temperature, although the effect is also not very drastic (Fig. |j). The shapes of 
both distributions resemble that of a Poissonian, as it should be in a grand canonical 
ensemble [24], but this Poissonian is slightly modified by resonance decays and isospin 



conservation. Compared with the empirical distribution of Fig. [j], all the curves of Figs. 
[3] and |] are much too broad. This means that the constraint on the partition function 
originating from isospin conservation alone is too weak, and that there are still too many 
open channels broadening the pion distribution. 

Now, by additional conservation of "external" quantum numbers a stronger constraint 
is put on the number of possible multi-meson channels. Due to the weak temperature 
dependence of the distribution and the large numerical effort, we have restricted the 
calculations to a single temperature value of T = 160 MeV. In order to account for the 
unknown ratio X r of the spin-state occupation of protonium just before annihilation, we 
have plotted the pion distribution for several values of W(0~ + ) in Fig. [|, using Bose 
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3 AVERAGE PION MULITIPLICITY AND MULTIPLICITY DISTRIB UTION 
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Figure 3: Pion multiplicity distribution with (n n ) = 5, for T = 100 MeV and T = 
200 MeV. The Boltzmann approximation is used and only the isospin and its third com- 
ponent are conserved. 
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Figure 4: Pion multiplicity distribution with (n n ) = 5, for T = 100 MeV and T = 
200 MeV with Bose statistics. Only isospin I and its third component I3 are conserved. 
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w(cr + ) = o.o 



(a) 



W (0" + ) = 0.2 



(b) 
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W(Q- + ) = <X8) 
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(e) 



W (0" + ) = 1.0 



(f) 
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Figure 5: Pion multiplicity distributions with (n n ) = 5, for T = 160 MeV using Bose 
statistics. In addition to the isospin conservation, also AMC and conservation ofP, C, and 
G parity are included. Six different values for the spin-state occupation of the protonium 
are taken into account, which are characterized by (a) W(0~ + ) = 0.0, (b) W(0 h ) = 0.2, 
(c) W(Q- + ) = 0.4, (d) W(Q- + ) = 0.6, (e) W(0~ + ) = 0.8, and (f) W(0~ + ) = 1.0. 
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4 CONCLUSIONS 



statistics. The behaviour of the calculated pion distribution with increasing W(0 + ) can 



numbered n n states dominate strongly, especially the n n = 4 and n n = 6 states. For 
W(0~ + ) = 0.8 (Fig. ||e) the states with = 3 up to n n = 7 are rather equally distributed. 
Finally, in the limit of exclusive occupation of the _+ state, the odd numbered states with 
n n = 3, 5, and 7 dominate the distribution, see Figure ||f. For statistical occupation of 
the spin states according to their combinatoric weight, W(0~ + ) = 0.25, we would have 
a result close to that of Fig. [|b. All in all, the obvious conclusion is that generally it 
is very difficult to get a Poissonian or even a Gaussian distribution when the external 
quantum numbers are conserved. This complete destruction of a smooth distribution 
by the conservation of external quantum numbers is clearly a consequence of the large 
fraction of directly produced pions in the partition function, which overwhelms the indirect 
production through resonance decays due to the larger phase space. We want to emphasize 
here that this is not an artefact of our grand canconical approach, but should also be 
observed in microcanonical calculations. 



4 Conclusions 

We have shown that, in the framework of a statistical model, Bose statistics and AMC 
together with conservation of the parity quantum numbers are crucial ingredients in the 
description of the multi- meson channels resulting from pp annihilation at rest. In such an 
approach, the projection operator formalism appears to be a convenient and consistent 
method in order to implement the conservation of non-Abelian quantum numbers. 

Isospin conservation alone yields a smooth but much too broad pion multiplicity dis- 
tribution, nearly independently of the temperature T, and the interaction volume V is 
unphysically big. Only Bose statistics together with conservation of external quantum 
numbers provide a reasonable interaction volume of V = [5 ... 20] fm 3 in the range of 
T = [140 . . . 200] MeV. The resulting pion distribution, however, shows no resemblance 
to a Gaussian or Poissonian: either the even or odd numbered states are strongly 
dominating. Thus, even in a statistical approach, once the effects of quantum number 
conservation are taken into account, the appearance of a smooth, Gaussian multiplicity 
distribution remains a mystery. That microscopic dynamical effects from the underly- 
ing QCD mechanism of annihilation should restore the apparent statistical nature of the 
empirical multiplicity distribution, is hard to believe. 



be summarized as follows. Up to a value 




even 
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In light of this dilemma, we suggest that the solution might be buried in the inclusive 
pion data themselves: the available bubble chamber data are rather old, and several ex- 
periments with probably different target densities have been exploited. For this reason, it 
would be desirable to make a new measurement of the pion multiplicity distribution at 
various fixed values of the hydrogen density under stable conditions. 
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